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Abstract 

We study an effective theory of QCD at high density in detail, including the 
finite temperature effects and the leading order correction in \f\x expansion. 
We investigate the Cooper pair gap equation and find that the color-flavor 
locking phase is energetically preferred at high density. We also find the 
color-superconducting phase transition occurs in dense quark matter when 
the chemical potential is larger than 250 ± 100 MeV and the temperature is 
lower than 0.57 times the Cooper pair gap in the leading order in the hard- 
dense-loop approximation. The quark-neutrino four-Fermi coupling and the 
quark- axion coupling receive significant corrections in dense quark matter. 

PACS: 12.20.Ds, 12.38.Mh, ll.10.Gh 

Keywords: Finite density/temperature QCD; Effective Theory; Cooper pair 



Typeset using REVTeX 



*On leave from Department of Physics, Pusan National University, Pusan 609-735, Korea. 
E-mail address: dkhong@hyowon.cc. pusan. ac.kr 



1 



I. INTRODUCTION 



Effective field theory is indispensable to describe physical processes at energy lower than 
the characteristic scale of theory, especially when the microscopic physics is too complicated 
to be used as in the case of chiral perturbation theory Jl| or just unknown as in the standard 
model. Since effective field theory has been very useful and successful in explaining the low 
energy data, it has now become an everyday language in physics . 

Color superconductivity in cold and dense quark matter || has been studied quite in- 
tensively in recent years by using effective field theory methods |5|-|§. Most important 
ingredient in those studies was the existence of effective four-quark interactions that mimic 
the four-fermion interaction generated by phonon exchange in BCS superconductivity, which 
are assumed to be generated in the effective Lagrangian either by exchange of massive glu- 
ons || or by strongly coupled instantons For low density quark matter, where strong 

interaction is no longer weak, such four-Fermi interactions may describe QCD at low density 
as model interactions, thus leading to color superconducting gap of QCD scale (~ Aqcd)- 
But, it was recently argued quite convincingly by Son that at least for high density quark 
matter, where perturbative QCD is applicable, magnetic gluons are not Debye-screened but 
only dynamically screened or Landau damped. Such long-range magnetic gluon interaction 
is more important and leads to a bigger Cooper-pair gap, A ~ fxg~ 5 exp(— 3n 2 /y/2g s ), than 
the usual BCS gap, which is subsequently confirmed by the Schwinger-Dyson analysis for 
the Cooper-pair gap 

Recently, an effective field theory of QCD at high density [|11|] is derived systematically 
in the power expansion of by integrating out the anti-quarks which are decoupled at low 
energy in dense quark matter. In the effective theory, four-quark operators are generated 
at the one-loop matching and the electric gluons are screened due to quarks in the Fermi 
sea while the magnetic gluons are not at least in perturbation |iO|Ji5H l8|. At scales below 
the electric screening mass, the relevant interactions for quarks are the coupling with mag- 
netic gluons and the four-quark interaction with opposite momenta. Both interactions for 
quarks in the color anti-triplet channel with opposite momenta are shown to lead to color 
superconductivity HlQjJllfl . 

According to the quark-hadron phase continuity in QCD, conjectured recently by Schafer 
and Wilczek ||19|| , the confining phase of nuclear matter at low density is complementary to 
the Higgs phase quark matter at high density. We support this conjecture by showing that in 
the effective theory, which is equivalent to QCD in the asymptotic density, the color-flavor 
locking diquark condensate is energetically more preferred. Furthermore, in the effective 
theory, the higher order corrections are systematically calculable and one can easily estimate 
quantities like the critical density, which are determined by the sub-leading operators in 1/fi 
expansion. 



In this article, we study the effective lagrangian of QCD at high density derived in [11 



in detail, including the higher order corrections in expansion and the finite temperature 
effects. We then analyze the Cooper-pair gap equation more rigorously and calculate the 
gap, the critical temperature, and the critical density. As applications, we examine the 
effect of marginal four-quark operators on the neutrino-quark four-Fermi coupling and the 
quark-axion coupling in super-dense quark matter. 
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II. HIGH DENSITY EFFECTIVE THEORY 



The idea of effective field theory is quite simple and its construction can be performed 
in a systematic way. First, one identifies or guesses the right degrees of freedom to describe 
the low energy processes. Then, one integrates out the irrelevant degrees of freedom, which 
results in nonlocal interactions. By expanding the nonlocal interactions in powers of mo- 
mentum at low energy, one derives a (Wilsonian) effective Lagrangian, which is usually done 
in practice by matching all the one light-particle irreducible amplitudes in the full theory 
with those amplitudes in the effective theory in loop expansion. Since the higher dimen- 
sional operators have smaller effects at low energy, the effective theories provide a very useful 
approach to describe low energy physics. 

A system of degenerate quarks with a fixed baryon number is described by the QCD 
Lagrangian density with a chemical potential //, 

£ QCD = i>%Pi> - \F%F a ^ + ^70^, (2.1) 

where the co variant derivative _D M = + ig s A a ^T a and we neglect the mass of quarks for 
simplicity. The chemical potential is introduced as a Lagrangian multiplier to constrain the 
system to have a fixed baryon number, N B . At zero temperature, all the states up to the 
Fermi surface are filled; 

where V is the volume, Nf the number of quark flavors, and f)p is the Fermi momentum. 
The chemical potential /i, which defines the energy at the Fermi surface, provides another 
scale for QCD. Here we consider a super-dense quark matter where quarks are very dense 
so that the average inter-quark distance is much smaller than the characteristic length of 
QCD, 1//J, <C 1/Aqcd- For high chemical potential, the pair creation of a particle and an 
anti-particle is suppressed at low energy because of the energy gap n, provided by the Fermi 
sea. Therefore, the relevant degrees of freedom at low energy will be just particles and holes, 
excited near the Fermi surface, together with gluons. 

We will derive a low energy effective Lagrangian for cold and dense quark matter by 
integrating out the high frequency or irrelevant modes, d la Wilson, as we run QCD from 
high energy E > fi to low energy E ~ Aqcd • Since QCD is asymptotically free, quarks will 
interact weakly at energy E > jj, (^> Aq CD ) and their spectrum will be approximately given 
by the free Hamiltonian: 

(d-p-ii)ip± = E±ip±, (2.3) 

where a = 707 and ip± denote the energy eigenfunctions with eigenvalues E± = —fi ± \p\, 
satisfying a ■ pip± = ± \p\ ip±, respectively. At low energy E < fj,, the states ip + near the 
Fermi surface, \p\ ~ fi, are easily excited, while ip-, which correspond to the states in the 
Dirac sea, are completely decoupled due to the presence of the energy gap fi provided by the 
Fermi sea. Therefore the right degrees of freedom below /i consist of gluons and *0 + only. 

If fi = 0, integrating out the high frequency modes in QCD just renormalizes the coupling 
constants without generating any new interactions as it should be for a renormalizable 
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theory like QCD. However, when fi 7^ 0, the energy spectra of particles and anti-particles 
are asymmetric. Anti-particles which correspond to the absence of states (iff-) in the 

Dirac sea, have energy larger than \x and will be decoupled at energy scale lower than fx. 
Therefore, the minial quark-gluon coupling that connects ?/>_ with ip + (or vice versa), has to 
be removed at low energy, generating a new interaction among ip + . We will see later that 
integrating out the modes whose frequency is higher than fi results in an effective Lagrangian 
which is drastically different from QCD. 

We now integrate out the antiparticles (V>_) or the states (ip~) in the Dirac sea to derive 
the low energy effective Lagrangian. We first decompose the quark momentum into the 
Fermi momentum and a residual momentum as 

Vy, = /«V + ln> ( 2 - 4 ) 

where the residual momentum |/^| < fi and = (0,v F ) with Fermi velocity v F . The 

magnitude of Fermi velocity, \v F \ = p F / \Jp 2 F + rn 2 , will be taken to be unity, for we assume 
the quark mass m — 0. 

Since we will integrate out the high frequency modes such that we will be left with 
excitations only near the Fermi surface at low energy, it is convenient to decompose the 
quark field by the Fermi velocity Q 

iP( x ) = e^^ivp, x), (2.5) 

Vp 

where 

WF,x)=f jPu^vp,^-^ (2.6) 

carries the residual momentum l^. We introduce projection operators P± = (1 ± a • v F )/2 
and define 

ip±(v F , x) = P±ip(v F , x). (2.7) 

In the limit 1/ fx — > 0, if)-(vp, x) corresponds to the states in the Dirac sea and ip + (vp, x) to 
the states above the Dirac sea. 

In terms of the new fields, the QCD Lagrangian for quarks at high density becomes 

£quark = $(%) (« P + ^7°) $(x) 

= ]T [$+ @F, x)ij°V ■ (vp, x) + ^_ (v F , x) 7 ° (2/i + iV-D) tfj_ (v F , x) (2.8) 

Vp 

+7p-(v F ,x)ij t [D fl ip + (v F ,x) + 4> + (v F ,x)i'j t tD^_(v F ,x) , 



1 For a given quark momentum, the corresponding Fermi velocity is determined up to 
reparametrization; v F — ► v F + 6l±/fi and I —* I — 51, where 5l± is a residual momentum per- 
pendicular to the Fermi velocity. As in the heavy quark effective theory [20|, the renormalization 
of higher-order operators are restricted due to this reparametrization invariance. 
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where = (l,v F ), V 11 = (1, —v F ), J± = — Ty 1 with '-y = (7 , vfVf ■ 7), and we have used 

iP + (v F ,x)^iP + (v F ,x) = ^ + P_ 7 M P+^ + = V^ +1 °iP + 

4>^ (v F , x)Yip- {vf, x) = ^_P + 7 M P_^- = V- (2.9) 

1p + (vF,x)'y fl 1p-(vF,x) = ljj + P^J^P^tp- = '0 + 7j^'0_ 

i)^{v F ,x)^ip + {v F ,x) = V>_P + 7 M P_^+ = V>_7^+. 



In Eq. ( |2.8| ) the incoming quark and the outgoing quark have the same Fermi velocity vp, 
since, at low energy, the momentum carried away by gluons can be compensated by the 
residual momentum of the quarks to conserve the momentum, without changing the Fermi 
velocity. 

As we can see from Eq. ( |2.8|) , the propagation of quarks at high density is 1+1 dimensional 
if l\i ' \i is negligible. This dimensional reduction at high density can be also seen if we 
decompose the quark propagator as following: 

iS F (p) 



[1 + ie)p°7° — p ■ 7 + /i7° 



= ~ ' ■, + | 9 —, V ■ • / + °<?Jti> ( 2 - 10 ) 
2 I -V + iel 2 2\i + I ■ V + iel 

where quark momentum p M = fiv^+l^ and the it prescription is chosen such that not only the 
anti-particles but also the holes correspond to the negative energy states moving backward 
in time. Note also that, if l\j \i — > 0, 1/2 ^7° = 1/2(1 +a-Vp) is the projection operator that 
projects out particles and holes while 1/2 ^7° = 1/2(1 — a-v F ) projects out the states 
in the Dirac sea (ip-)- Since the excitation of quarks along the perpendicular direction to the 
quark velocity does not cost any energy in the leading order, the perpendicular momentum 
is nothing but an index that labels the degeneracy in quark just as the degeneracy in the 
Landau level under external constant magnetic field is labeled by the electron momentum 
perpendicular to the external magnetic field. 

At tree-level, integrating out the irrelevant modes ip_(vp, x) is tantamount to eliminating 
tp^(v F ,x) by the equations of motion, given as 

• • 00 / • \ n 



where D\\ = V^D^ and p±_ = ^D^. Plugging Eq. ( 2.11 ) into the Lagrangian for quarks 



Eq. ( |2.8| ), we obtain the tree- level effective Lagrangian of QCD at high density 

i/j + (vf, x)iV ■ Di/) + (v F , x) — ijr + 



4 

Vp 



1 

n 




III. ONE-LOOP MATCHING AND A FOUR-FERMI OPERATOR 

As we have seen in the previous section, integrating out the irrelevant modes generates 
a new coupling between quarks and gluons. This can be also seen diagramatically as shown 
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in Fig. 1: In the leading order, the propagator of fast modes ip- is replaced by a constant 
matrix 

and the exchange of ?/>_ generates a new interaction given as, using P_7^P_ = P_7j^, 

4°ff 3 -7T E i>+ A7o + " " " , (3.2) 

where the ellipsis denotes terms containing more powers of gluons and derivatives. By writing 
the interaction in a gauge-invariant fashion, we recover the interaction terms in Eq. ( |2.12| ). 



The effect of integrating out fast modes not only generates new interactions at low 
energy but also gives rise to quantum corrections to the tree-level couplings which can be 
calculated explicitly by matching the loop amplitudes. Usually the new interactions arise at 
tree-level upon integrating out the fast modes, but sometimes they may arise at quantum 
level, especially when they are marginal, as in the low energy effective theory of QED under 



external (strong) magnetic field |fJT 



As fermions under external fields f2T^p3f, in the presence of a Fermi sea || the scaling 



dimension of fields at low energy changes due to the change in the particle spectrum; the 
scaling dimension of quarks becomes —1/2 instead of —3/2, the canonical scaling dimension 
of fermions in 3+1 dimensions. The scaling dimension of fields is determined by the kinetic 
term in the action, which is for quarks in dense quark matter given in the high density limit 

as 

So = E/ tJ4i \f+$F, (h - Ml) M*f> + 4>-@f, I) (2fi + l- V) ^-{v F) l)] , (3.3) 



VF 



where l\\ = Vp • I is the residual momentum parallel to the Fermi velocity. Under a scale 
transformation Iq — > slo and L — > slu with s < 1, the quark fields transform as ip + (vF, I) — > 
s~ 3 / 2 ip + (vF, I) and ip^(vp, I) — > s~ 1 ip_(vp, I)- As we scale the momentum toward the Fermi 
surface or we decrease s — > 0, tp + scales like free massless fermions in (l+l)-dimensions, 
while ?/>_ scales like heavy fermions in (l+l)-dimensions and becomes irrelevant at low 
energy. Since four-Fermi interactions are marginal for (l+l)-dimensional (light) fermions, 
we expect that a marginal four-quark operator will arise in the high density effective theory 
of QCD by quantum effects if it is absent at the tree-level. To see this, we consider one-loop 
matching of a four-quark amplitude for scattering of quarks with opposite Fermi velocities, 
shown in Fig. 2, because only when the incoming quarks have opposite Fermi momenta the 
four-quark operators are marginal as we scale toward the Fermi surface 0. Since the effective 
theory amplitude is ultra-violet (UV) divergent while the amplitude in full QCD is UV finite 
by the power counting, we need a UV counter term for the one-loop matching in the effective 
theory, which is nothing but a four-quark operator. The one-loop four-quark amplitude of 
our interest in the effective theory is, reverting the notation if) for ip + henceforth, 

(#3, ^3; U4, k\ i> A±l° A±^(vf, x) ■ tfj /4_l7° / 4_ l V(^ ) y) |#i> h] v 2 , h) 

Jx,y 
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4/i- 



■(27r) 4 5 4 (^^)^(p3)7^7 7:T a TV(pi) ■ Hp 4 hil°llT c T d ^(p 2 ) 



(3.4) 



x 



(5 ac 5 b % p g ua + 8 ad 5 bc gfi{T g vp 



—i 



16tt 2 



1 , (-1 2 -16 



where i7j's are the Fermi velocities of incoming and outgoing particles with vi = —V2, p% = 
fiVi + k, 7 the Euler number, Yf Pi = P3+P4 — Pi~ P2, I = h~h, and A is the renormalization 
point. 

In QCD, the quark-quark scattering amplitude is 



^ CD = |(P3,P4| 



5^ / ip(v F ,x) fi ± (x)ip(v F , 

Jx 



X) 



\Pl,P2) 



gt(2n)H\J2'P^(P*h±rt TaTb ^P^ (3-5) 



where 



W 

Jo 



1 



1 



, 2fM + V ■ (Is + q) 2fi + V ■ {h ~ q) (q + l) 2 q 2 
1 1 1 

, 2/i + V • (/1 - q) 2/i + V ■ (k -q){q + l) 2 q 2 



(3.6) 
(3.7) 



In the limit that the residual momenta U — > 0, we find I\ = I2, if we rotate the qo axis into 
v F ■ q axis, and thus the color and Lorentz structures of the amplitudes in both theories are 
same. To perform the integration, we take, for convenience, the external residual momenta 
are perpendicular to the Fermi velocity, V • 1% = 0. Then, we get for lij \x — > 



h 



dx 



l-x 



64n 2 fi 2 



dy 
2 - In 



[(l-x- y)q 2 + x(q + I) 2 + yq 2 - y(v F ■ q + 2/i) 1 
-I 2 ' 



(3.8) 



4/x 2 



Note that the scattering amplitude A® CD is IR divergent as U — > 0, which is the same infrared 
divergence of the amplitude in the effective theory, Eq. ( |3.4| ), as it should be because the 
low energy effective theory is equivalent to the microscopic theory in the infrared limit by 
construction. Since both amplitudes have to be same at the matching scale A = /i, we need 
a four-quark operator in the effective theory, given by the difference of the amplitudes at 
the matching scale [l: 



.QCD A eS 



64ttV 



1 



- - 7 - 4 (27r) 4 5 4 (£ ' K )^( ft )7KrTV(pi 



x^fo) (l,luT a T b + lvl ,T b T a ) V»(p 2 ) + 0(k/v) 
(P3,Pi\ i J d 4 x C\ f \pi,p 2 ) + 0(li/fi), 



(3.9) 
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where, at the matching scale /i, the new effective four-quark operator generated at one-loop 
is given as 



C 1 



4/ 



„bare 

yi 

' 2/i 2 



x ^(-ifr, x) (lfj,luT a T b + lvlll T b T a ) V(- 



(3.10) 



In the modified minimal subtraction, the renormalized four-quark coupling is gj en = 2a 2 
at the matching scale \x. It is convenient to rewrite the marginal four-quark operator by 
Fierz-transforming the product of gamma matrices and SU(3)c generators in the amplitude 
Eq. ( ft.lOQ . Using {T a ) tu (T a ) vs = l/2S ts S uv — l/65t u S vs and the Fierz transformation, we get 



[P+llllP+h { TaTb ) ut ■ [ P - (l^±T a T b + l v i_l,i_T h T a ) P_ 

"1 



lm;vs 



[P+lll±P+]ji (T a T b ) 



at 



nab 



d abc T c ) + {P_a (U ,PJ) f abc T c 



Im.vs 



Int. 
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6 



SutSvs + -5 US S t v ) + (P+l5)ji (-P-75) 



lira 



— 8 u t5 vs + — S us 5t v 



(3.11) 



(P+)ji (P-)lm ^ g ^uv;ts 



13 
9\71 



5 s 

w uv:ts 



(P+75)^ (P-To) lm (V25t,u ~ ^uv^ 



where i,j,l,m denote the Dirac indices and t,s,u,v the color indices. In the last line, 
the color indices in the operator are further arranged to decompose the amplitude into 
the irreducible representations of SU(3)c by introducing invariant tensors in color space, 
$uv,ts = (Ktfivs + KsSvt) / V2 and 5£ v . t3 = (5 ut 5 vs - 5 us 5 vt ) / y/2. 
The four-quark operator in the effective theory becomes then 



r l - 



\l H [9us;tv1pt(VF, x)lp s (v F , X)lpl(-V F , x)lp u (-V F , x) 
v F 

+h us; t v ipt(vF, x)j 5 ip s (v F , x)ipl(-v F , x)^ u {-v F , x) 



(3.12) 



with g us . tv = g^ s . tv - g 6 5^ s . tv , h us;tv = h^ s . tv - h 6 6% s s . tv . The value of couplings at the 
matching scale fi are given as g% = 4V2a 2 s /9 = 20 6 /(13) and = A^2a 2 s = 2h 6 . 



IV. SCREENING MASS 



As discussed in [16-18,24 



the quark loop correction to the vacuum polarization tensor 
gives rise to color screening in quark-gluon plasma, while the gluon loop renormalizes the 
color charge. We first calculate the vacuum polarization tensor in QCD and then compare 



2 If we had matched a quark-quark scattering amplitude with momenta not opposite to each other, 
we would get similarly a four-quark operator in the effective Lagrangian. But, since it is irrelevant 
as we scale toward the Fermi sea, it will not affect the low energy dynamics significantly. 
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it with the one in the effective theory when we perform the one-loop matching for the gluon 
two-point amplitude. 

Since the Feynman propagator for quarks of mass m in matter at zero temperature is 
given by 



iSp(x) 



d A p 



-ip-x 



(27r) 4 (1 + ze)po7° ~ P ' 7 + ^7° ~~ m 
we find, performing the integration over po, 

d 3 p $ + 171 



(4.1) 



iSp(x) = 9(x 
- 0(-x 



(2tt)3 2p 
d 3 p 



-6(p - fi)e 



-ip-x+ifixo 



(27T) 



+ m 



2p 



9(p-p )e 



-ip-x+i^iXQ 



+ 



— m 



tp-X + ifjLXQ 



2p 



(4.2) 



where Po = y\p\ + m 2 . We see that this propagator agrees with the one obtained by the 
canonical quantization |T7| except the overall phase factor due to the shift in the energy by 
p to set the energy of the Fermi surface to be zero. 

In terms of this full propagator Eq. (3~I), the quark-loop contribution to the one- loop 
vacuum polarization tensor becomes 



£Jd A xe-*'{J£(x)Jt;(0)) 
d 3 q 1 1 



iN f 2, 



(2vr) 3 2q 2k 



(4.3) 



qo + h+Po- it 



T^(q, k! 



qo-h + Po- ^ 
0(k - y) 



qo - h + Po + it 



-qo - h + po + ie 



T^ u (q, k") 



(4.4) 



where T M ^(g, k) = Tr py M <^ r y v ft], k = q + p, k' = — q — p, and k" = —q + p. 

If we rewrite the step functions as 9(qo — p) = 1 — 9(p — qo) and 9(k — p) = 1 — 9(p — ko) 
for the last two terms in Eq. ( |4.4| ), one can easily see that the quark-loop contribution of the 
vacuum polarization tensor consists of two parts, one due to the matter and the other due 
to the vacuum, FQ ull (p) = U^ m£Lt (p) + U^ vac (p), where Il^ac is the quark- loop contribution 
when there is no matter (i.e. p = 0), 



abvac 



tN f 2, 



d 3 q 1 1 
(2^F2^2^ 



T^(q, k') 



T» u (q,k") 



qo + ko+Po- it -qo - k Q + p + it 



(4.5) 



where Nf is the number of light quark flavors. For <C p, the matter part of the vacuum 
polarization becomes, with M 2 = Nfg 2 p 2 /(2ir 2 ), 



afemat 



iM 2 



-5, 



ab 



dQ 



-2p-v F V>*V u uv WV + WV"^ 
+ 9^ n 



An V p ■ V + iep ■ Vp 



(4.6) 



which is transversal, P^abmatip) = EI 



Now, we calculate the quark contribution to the one-loop vacuum polarization tensor in 
the effective theory. Quarks in the (high density) effective theory are almost on-shell at low 
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energy and the quark current, which consists of states near the Fermi surface, is proportional 
to its velocity, V 11 = (l,Vp): 



J^(x) = J2Mv F ,x)^T a MvF,x) = J2V^ + (v F ,x) 1 °T a MvF,x), (4.7) 



Vp Vp 



where we used P + tp + (vp,x) = ip + (vp,x) and P_7 M P + = P_7 ( V M . Therefore, the gluons 
transversal to the quark velocity is decoupled at low energy and the gluons longitudinal to 
the quark velocity is coupled to the quark color charge density in a combination given as 

Al{x)J^{x) =J2MvF,xh°T a ^ + (v F ,x) V-A a (x). (4.8) 

Vp 

Therefore, in the leading order in expansion, the quark-gluon coupling does not involve 
the spin of quarks like a scalar coupling. For a given quark, moving nearly at a Fermi 
momentum fiVF, only one component of gluons, namely Aq — vp ■ A combination, couples to 
the quark in the leading order. 

With the quark current given in Eq. Q4.7p , we obtain the quark-loop contribution to the 
vacuum polarization in the effective theory as 

KM = 9s J d 4 x e^^(0\TJW F ,x)Ji;(# F 0)) (4.9) 

Vp 



Vp 



a - Vp \ n % r 

7 / , • 7 



-— s ab y v^v 1 - 

8tt 4^ 

Vp 



1 - 



(p° +p- Vp) 2 



M , (4.10) 



where = (p° , vpvp ■ p) . We note that the quark- loop contribution to the vacuum polariza- 
tion tensor in the effective theory is exactly same as the sum of the first and second terms 
in the QCD vacuum polarization, Eq. ( |4.4j ). This is what it should be, since the effective 
theory does reproduce the contribution by quarks and holes in QCD when the external gluon 
momentum p — > 0. The third and fourth terms in Eq. (|4.4|) , due to the quark and anti-quark 
pair creation, is absent in the effective theory because the anti-quarks are integrated out. 
The effect of anti-quarks will be taken into account when we match the gluon two-point 
amplitudes. 

To match the gluon two-point amplitudes in both theories, we therefore need to add a 
term in the one-loop effective Lagrangian, 



M 2 

(4-11) 



Vp 



which also ensures the gauge invariance of the effective Lagrangian at one-loop. The static 
screening mass can be read off from the vacuum polarization tensor in the limit po — ► 0, 
which is in the effective theory 

KUPo - 0,P) - -iM 2 6 ab 6^°. (4.12) 
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We see that the electric gluons, Aq, have a screening mass, M, but the static magnetic 
gluons are not screened at one- loop due to the added term Eq. ( 4.11| ), which holds at all 
orders in perturbation as in the finite temperature p^JT5[| . 

Finally, by matching the quark two-point amplitudes, we get a one-loop low-energy 
(Wilsonian) effective Lagrangian density; 



1 2 A/f 2 - 1 

C cS = --(1 + a x ) (F; v ) - j^ATA^ + (1 + h)^D^ - — (1 + Cl )^ ( 7 . 



2(i 



2fi 2 



(finite - 9e€s;tv) *Pt(vF, x)tfj s (v F , x)ipl(-v F , x)ip u (-v F , x) 

+ ( h ^us;tv - he S us;tv) ^4 (#F, ^)l^s (v F , x) (~V F , x)j 5 <lp u (~V F , x) 



(4.13) 



+ ■ ■ 



where the summation over v F is suppressed and the coefficients ai,bi,c\ are dimensionless 
and of order The ellipsis denotes the irrelevant four-quark operators and terms with 

more external fields and derivatives. 



V. RG ANALYSIS AND GAP EQUATION 

As we scale further down, the effective four-quark operators will evolve together with 
other operators, which can be seen by further integrating out the high frequency modes, 
Sfi < < ft. The scale dependence of the four-quark operators has three pieces. One is 
from the one-loop matching condition for the four-quark amplitudes and the other two are 
from the loop corrections to the four-quark operators. Putting all contribution together, we 



find the one-loop renormalization group equations for the four-quark operators |TTJ to be 



where % = (6,3), ge = —g^, g^ = g% and % = (a/2/9)(13/4, 1/2) and 5i = (—1,2). Since 
in the high density limit the quark-gluon coupling is (l+l)-dimensional, the quarks do not 
contribute to the running of the strong coupling. The one-loop (3 function for the strong 
coupling constant at high density is (3{a s ) = — ll/(27r)a: 2 . 

To solve the RG equation, Eq. ( |5.1| ), we introduce a new variable 

1 /_ 7rln2 e , A .\ , 11 . . . 

yi =22^)[ 9 ' i + -r 6iasiA) ) and t = ^ a ^ s (5-2) 

Then, the RG equations becomes 

t - - jr+w (5 - 3) 

where a< = (1325; In 2 + 144^ - 5 2 )/(17424). Further, letting y { = /i(t)/(l + t), we get 



Integrating Eq. (|5.4j ), we get 
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f m _ ~ + tr {Cl - C2) - gi(/,(o) - c 2 ) 

M) mo)-ci)(i+t)-(^-<^-ifi(o)-c 2 ) ' 1 J 



where Ci = (1 + VI - 4a,)/2 = 1 - C 2 and £(0) = (&(//) + tt In 2<5 J a s (/z)/6) /(22ttq! s (//)). 
Since a; < 1, we take C x ~ 1, C 2 < 1. And also, for /i > A QC d, /i(0) ~ 5* In 2/(132) < 1. 
Therefore, for 1 + t — > + or s — > exp [— 27r/(llot a (/i))], we get ^(t) + 7rln2<5jQ; s (t)/6 ~ 
227rC 2 a s (t). Namely, 



&(A) ~ ^a a (A) 



(ln2) 2 2 

71 - ~uT 5i 



(5.6) 



At a scale much less than the chemical potential, A <C fi, g&(A) ~ 0.29a s (A) and (73(A) ~ 
0.04a s (A). Similarly for hi, % = ( v / 2/2)(2, 1) and $ = (-1, 2) and we get h 6 (A) = 0.81a fl (A) 
and /13(A) ~ 0.4a s (A). 

At a scale below the screening mass, we further integrate out the electric gluons, which 
will generate four-quark interactions. For quarks moving with opposite Fermi momenta, the 
electric-gluon exchange four-quark interaction is given as 

C l9 3 - 9 ^^Y.^° T ^F,x)i, l0 T a ^-v F ,x). (5.7) 

Vp 

Using Tt u T£ s = l/28 ts S uv — l/65 tu S vs , we find that the four-quark couplings are shifted 
as g- 3 (M) - g- 3 (M) + 2^ s 2 (M)/3 ~ 0.95£ 2 (M) and g 6 (M) - g 6 (M) + \/2<7 2 (M)/3 ~ 
0.49g 2 (M), while h^s are unchanged since Eq. (|5.7|) does not involve 75. 

As we approach further to the Fermi surface, closer than the screening mass M, the four- 
quark operators in the color anti-triplet channel become stronger, because the (3 function 
for the attractive four-quark operators is negative, (3{g 3 ) = — g|/(47r 2 ). If the four-quark 
interaction is dominant at low energy, it leads to vacuum instability in the infrared region 
by forming a color anti-triplet condensate or Cooper pair, But, since the long-range color- 
magnetic interactions also become strong at low energy, we need to consider both interactions 
to determine the Cooper-pair gap. 

Since both relevant interactions are attractive for a pair of quarks in color anti-triplet 
channel with opposite Fermi momenta, they may lead to condensates of quark pairs in 
color anti-triplet channel. To describe the Cooper-pair gap equation, we introduce a charge 
conjugate field, 

C0c)i {VF, X) = Cijtpjt-VF, x), (5.8) 

where i and j are Dirac indices and the matrix C satisfies C~ 1 7 At C = — 7J. Then, we can 
write the inverse propagator for ty(v F ,x) = {ip{vp,x) ) ip c {vF ) x)) T as 

S ( ^' /) = 7o l -At(/„) Z( h )i v)> (5 ' 9) 

where Z(h) is the wave function renormalization constant and A(Zy) is the Cooper-pair gap. 
The Cooper-pair gap is nothing but a Majorana mass for quarks. When Cooper-pairs form 
and Bose-condense, the quarks get a Majorana mass dynamically and the Fermi sea opens 
up a gap. 
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Schwinger-Dyson (SD) equations are infinitely coupled integral equations for Green func- 
tions. In order to solve the Schwinger-Dyson equations, we need to truncate them consis- 
tently. Since g s (M) is still weak at high density M ^> Aq CD , one may use the ladder 
approximation which is basically the weak coupling expansion of the SD equations, consis- 
tently with the Ward-Takahashi identity. Another gauge- invar inat truncation is so-called 
hard-dense-loop (HDL) re-summation, which can be derived gauge invariantly from the 
transport equation |T7] as in the finite temperature Po| . When the gluon momentum is of 



order of the screening mass M ~ g s fi, the quark one-loop vacuum polarization is also same 
order as we can see from Eq. ( 4.1(Jp . Therefore, if the important momentum scale of the 



diagrams is of order of M, one needs to re-sum all the bubble diagrams to get the consistent 
gluon propagator. Since the HDL re-summed gluon propagator correctly incorporates the 
medium effects like screening at high density, we will solve the SD equations for the quark 
propagator in the HDL approximation to calculate the Cooper-pair gap. 

In the HDL approximation, the SD equations for the quark propagator are given as, 
neglecting small hi four-Fermi couplings, 

zin) - 1] , • v = H*>" fv»o^ - ^W A gT + * (5.io) 

A(P„) = (-^/^(p-^^^ + l/^^, (5.1!) 

where the gluon propagator is given in the HDL approximation as,| following the 
notations used by Schafer and Wilczek 



where £ is the gauge parameter and the projectors are defined by 



PS = ^~ k A ^o = = PoT (5-13) 



1*1 



P^ = -9»» + ^-P^ (5.14) 

Before solving the SD equations, we first show that the preferred solution of the SD 
equations supports the color-flavor locking condensate predicted in 0. Since the gluon 
interaction is vectorial, the gluon exchange interaction in the gap equation does not dis- 
tinguish the handedness of quarks and thus it will generate same condensates regardless of 
handedness; \{4)l'4 , l)\ = \{'^r'4 , r)\ = KV'iV'i?)!) suppressing other quantum numbers. But, 
the four- Fermi interaction in the effective Lagrangian, Eq. (|4.13|) , is chiral, which can be 
seen if we rewrite it in the handedness basis as following; 



3 In the Schwinger-Dyson equation the loop momentum should take the whole range up to the 
ultraviolet cutoff, which is the chemical potential /i in the case of high density effective theory. 
Hence the gluon propagator includes both magnetic and electric gluons. 
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r l 93 + h rA 

■*"4f ^ o. 2 °us;tv 



2/j? 

+ 93 - h 5 A 



us;tv 



ilit^F, x)ip Ls (v F , x)ipl v (-v F , x)tjj Lu (-v F , x) + (L^ R)\ (5.15) 
tp{ t (v F , x)ip Ls (v F , x)^ Rv {-v Fl x)ip Ru (-v F , x) + (L^ R) . 



We see that the LL (or RR) four- Fermi coupling is bigger than the LR four- Fermi coupling, 
since > 0. Therefore, the gap in LL or RR channel will be bigger than the one in LR 
channel due to the difference in the four-Fermi couplings. Thus, the LL or RR condensate 
is energetically more preferred than the LR condensate. We also note that in the effective 
theory the gluons are blind not only to flavors but also to the Dirac indices of quarks 
because they couple to quark currents in the combination of V ■ A like a scalar field, as given 
in Eq. (|4.8| ). Therefore, in the SD equations, the diquark Cooper-pair can be decomposed 
into color anti-triplet and color sextet. But, since quarks of opposite momenta are attractive 
in the anti-triplet channel while repulsive in the sextet channel, without loss of generality, 
we can write the Cooper-pair gap in color anti-triplet Since quarks are anti-commuting, 
the only possible way to form diquark (S-wave) condensate is either in spin-singlet or in 
spin-triplet: 



^Li a (v F} x)lp L b jf3 (-V F} x) 



lf>Ria(v F , x)lp R b jp (-V F , X) 
abcj^ 



(Pi 



(5.16) 
(5.17) 



where a,b,c = 1,2,3 are color indices, a, /3, 7 = u,d, s, ■ ■ ■ , Nf flavor indices, and i,j = 
1,2 spinor indices. Indices in the bracket and in the curled bracket are anti-symmetrized 
and symmetrized, respectively. But, the spin-one component of the gap, K{ a/3 } c , vanishes 
algebraically, since i[)(v F , x) = 1/2 (1 + a • v F ) i/j(v f , x) and (1 + a • v F )u(l — a ■ v F )ij = 0. 
When Nf = 3, the spin-zero component of the condensate becomes (flavor) anti-triplet, 



K[ a p] c {p F ) = e a/ 3 7 K2(p F ). 



(5.18) 



Using the global color and flavor symmetry, one can always diagonalize the spin-zero con- 
densate as K] = S^K^. To determine the parameters, K u , Kd, and K s , we need to minimize 



the vacuum energy for the condensate. By the Cornwall- Jackiw- To mboulis formalism |27 
the vacuum energy in the HDL approximation is given as 



V(A) = -Tr In S~ L + Tr In + Tr (S 



-1 



47T 



2 9 



E 



d 2 i\\ 
W) 



In 



11 + A?(J B 



+ 



(2PI diagrams) 
A? (ill) 



1 

2 "l\ + A?(Z|, 



+ h.o., 



(5.19) 



where h.o. are the higher order terms in the HDL approximation, containing more powers of 
coupling g s , and Aj's are the eigenvalues of color anti-symmetric and flavor anti-symmetric 
9x9 gap, A^. The 2PI diagrams are two-particle-irreducible vacuum diagrams. There is 
only one such diagram (see Fig. 3) in the leading order HDL approximation. 



4 At high but finite density, the Cooper-pair gap contains a small component of color-sextet [26| 
But we will ignore this for simplicity. 
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Since the gap depends only on energy in the leading order, one can easily perform the 
momentum integration in (|5.19|) to get|^, 



2 



^-0-43^El^(0)| 2 , (5.20) 
where in the second line we used an approximation that 

Were A» independent of each other, the global minimum should occur at Aj(0) = const, for 
alH = 1, • • • , 9. But, due to the global color and flavor symmetry, only three of them are 
independent. Similarly to the condensate, the gap can be also diagonalized by the color and 
flavor symmetry as 

Aff = e^A.SJ. (5.22) 

Without loss of generality, we can take |A U | > |A d | > |A S |. Let A d /A u = x and A s / A u = y. 
Then, the vacuum energy becomes 

V(A) ~ -0.43^ |A U | 2 f(x, y), (5.23) 

where f(x, y) is a complicate function of — 1 < x, y < 1 that has a maximum at x = 1 = y, 
f{x,y) < 13.4. Therefore, the vacuum energy has a global minimum when A u = A^ = A s , 
or in terms of the eigenvalues of the gap 

A; = A u ■ (1, 1, 1, -1, 1, -1, 1, -1, -2). (5.24) 

Among nine quarks, ij}", eight have (Majorana) mass A u and one has mass 2A U . 

Since the condensate is related to the off-diagonal component of the quark propagator 
at high momentum as, suppressing the color and flavor indices, 



{^{v F ,x)^{-v F ,x)) ~ lim f -^\-J Hx ~ y) - A(/||) 



v-*xj (2vr) 4 lj-A 2 {l ]{ ) 



lim 

y->x 



47r z xii — y\\ 7m 



(5.25) 



where 7 m is the anomalous dimension of the condensate and the ellipsis are less singular 
terms. Being proportional to the gap, the condensate is diagonalized in the basis where the 



5 If the condensate forms, the vacuum energy due to the gluons also depends on the gap due to 
the Meisner effect. But, it turns out to be subleading, compared to the quark vacuum energy; 
V g (A) ~ M 2 A 2 ln(A//i) ~ g s ^ 2 A 2 (28]. 
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gap is diagonalized. Thus, we have shown that in the HDL approximation the true ground 
state of QCD with three massless flavors at high density is the color-flavor locking phase, 
K 1 = K for all 7 = u,d, s. The condensate takes 



i^Liaivp, x)ip Ljp (-v F , x)j = - (ijRi a (vF, x)if; R0 (-v F , x)j = e afSI K(p F ), (5.26) 

breaking the color symmetry, U(l) em , the chiral symmetry, and the baryon number symme- 
try. The symmetry breaking pattern of the CFL phase is therefore 

SU(3) C x SU(3) L x SU(3) R x U(l) cm x U(1) B ^ SU(3) V x U(1)q x Z 2 , (5.27) 

where SU(3)v is the diagonal subgroup of three SU(3) groups and the generator of U(1)q 
is a linear combination of the color hypercharge and U(l) em generator, 



Q = cos 9Q cm + sin 9Y 8 , 



(5.28) 



where tan# = e/g s . 

Now, we analyze the SD gap equation to see if it admits a nontrivial solution. Since the 
color-flavor locking gap is preferred if it exists, we may write the gap as 



Aab _ c abl \ 

Then the gap equation becomes, neglecting small hi couplings, 



(5.29) 



A(P||) = {-ig s f 



d 4 / 



■D, u (p-l)V 



,T"A(/||)(T") 



a\T 



-V u + i 



9s 



(2ir) 4 ~^ ^-A 2 (/,|) > 2 7 (2tt) 4 / 2 - A 2 (/,|) 



/ 



d 4 / A(/,|) 



, (5.30) 



where we use the bare vertex and take Z{p») = 1 in the leading HDL approximation. 
In the weak coupling limit, \ko\ <C \k\ and thus 



/' (/«'()■ />') - M 2 , G(k ,k) ~ ^-M 2 % 



\k\ 



(5.31) 



Since the gap has to be fully antisymmetric in color indices, we get 
Tt u A uv (T a )1 s = (^-S tv S us — -S tu 8 v ^j A uv = —-A ts 
After Wick-rotating into Euclidean space, the gap equation becomes 



(5.32) 



d 4 q 
(2k)' 



~3 9s 



V-P T -V 



(p-q)l + Ql + jM 2 \po-qo\/\p-g] 

.(p-?)S' 



93_ 



A(g||) 
?! + A 2 (gn)' 



(5.33) 



Note that the main contribution to the integration comes from the loop momenta in the 
region q 2 ~ A 2 and \q±\ ~ M 2 / 3 A 1 / 3 . Therefore, we find that the leading contribution is by 
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the first term due to the Landau-damped magnetic gluons. For this momentum range, we 
can take \p — q\ ~ \q±\ and 



V ■ P 1 ■ V 



-1 + 



/ A 4 / 3 \ 



(5.34) 



We also note that the term due to the four-Fermi operator is negligible, since g% ~ g\ at the 
matching scale /i. 

Neglecting (p — g) 2 in the denominator, the gap equation becomes at the leading order 
in the weak coupling expansion and l//i expansion 



m) 3 J (2tt) 4 



+ 



1 



gl + fM 2 |p -go|/|gi| gl + M 2 



ip ~ g)f 



A(gn 



gf + A 2 (g,|)' 



(5.35) 



The gl integration can now be performed easily to get 



97t7 (2vr) 2 g 2 + A 2 



In 



fM 2 | Po -go| 



+ 2 ln [W + 2^ 



(5.36) 



We see that in this approximation A(p\\) ~ A(p ). Then, we can integrate over vp ■ q to get 



Afro) 



36tt 2 J- 



, A (g ) 
ago , = hi 

" V% + A 2 V bo -go I 



(5.37) 



where A = • (/i/M) 5 e 3/2? . If we take A ~ A(0) for a rough estimate of the gap, 

2 



1 



367T 2 



In 



A 



or A ~ A exp — 



67T 

9 s , 



(5.38) 



As was done by Son [|H]], one can convert the Schwinger-Dyson gap equation ( |5.37[ ) into a 
differential equation to take into account the energy dependence of the gap. Approximating 
the logarithm in the gap equation as 



In |p - go | 



ln|p | if |po| > M, 
In | go | otherwise, 



we get 



pA"(p) + A'{p) + 



2a s A(p) 



9tt vV + A 2 



0. 



(5.39) 



(5.40) 



with boundary conditions pA' = at p = A and A = at p = A, where p = p . When 
p < A(p), the equation becomes 



(5.41) 



where r 2 = 2g 2 /(97r 2 ) and |A| is the gap at p — 0. We find A(p) = \A\J [ryp/\A\J for 
p <C || A|. When p ^> A, the differential equation ( |5.40| ) becomes 
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r 2 A 
4 p 



(5.42) 



which has a power solution, A ~ p ±vr ^ 2 . Since the gap vanishes at p = A, we get for p ^> A 

(5.43) 



A(p) = fisin f~ln-j . 



By matching two solutions at the boundary p = |A| we get 

B ~ |A| and |A| = Ae^ /r . 
The gap is therefore given as at the leading order in the weak coupling expansion^ 

3tt 2 



IAI 




V2gJ ' 



(5.44) 



(5.45) 



where c = 2 7 7r 4 iVy 5//2 e 3 ^/ 2+1 . This agrees with the RG analysis done by Son |H| (see also p9| ) 
and also with the Schwinger- Dyson approach in full QCD [p^ -|T4]|, The l/g s behavior of the 
exponent of the gap at high density is due to the double logarithmic divergence in the gap 
equation ( 5.33j ), similarly to the case of chiral symmetry breaking under external magnetic 
fields p2||30| , pT| . In addition to the usual logarithmic divergence in the quark propagator as 
in the BCS superconductivity, there is another logarithmic divergence due to the long-range 
gluon exchange interaction, which occurs when the gluon loop momentum is colinear to the 
incoming quark momentum (q± — > 0). 



VI. TEMPERATURE EFFECTS AND HIGHER ORDER CORRECTIONS 

When the quark matter is not very dense and not very cold, the effects of finite temper- 
ature and density become important. In this section we calculate the critical density and 
temperature. First, we add the corrections to the gap equation Eq. ( |5.30|) to see how 
the formation of Cooper pair changes when the density decreases. As derived in [IT], the 
leading corrections to the quark-gluon interactions are 



2n 



4[i 



(6.1) 



In the leading order in the HDL approximation, the loop correction to the vertex is neglected 
and the quark-gluon vertex is shifted by the l//z correction as 



-ig s v F - ig s 



(6.2) 



6 The gauge-parameter dependent term is subleading in the gap equation ( 5.33| ). Since the gap 
has to be gauge-independent, the gauge parameter dependence in the prefactor will disappear if 
one includes the higher order corrections. 
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where U is the momentum carried away from quarks by gluons. Then, the gap equation 
( |5.35| ) becomes 



Hp\\ 



2gi 



dH 

(2tt)< 



+ 



|/l| 3 + (7r/4)M 2 |/ -po| l\ + M* 



^ 2 + , f(l-^ 2 ) 



I 4 



Ago) 
/ 2 + A 2 



. (6. 



For a constant gap approximation, A(pn) ~ A, the gap equation becomes in the leading 
order, as p — > 0, 




/^ + A 2 



(6.4) 



Therefore, we see that, when // < \i c ~ e 3 A, the gap due to the long-range color magnetic 
interaction disappears. Since the phase transition for color superconducting phase is believed 
to be of first order f32| , |3~3"|| , we may assume that the gap has the same dependence on 
the chemical potential fi as the leading order. Then, the critical density for the color 
superconducting phase transition is given by 



H c = e n c exp 



4tt\/3 



(6.5) 



Therefore, if the strong interaction coupling is too strong at the scale of the chemical po- 
tential, the gap does not form. To form the Cooper pair gap, the strong coupling at the 
scale of the chemical potential has to be smaller than <7 s (/i c ) = 7r 2 /\/2. By using the one- 
loop P function for three light flavors, (3(g s ) = — 9/(167r 2 )g 3 , and the experimental value 
for the strong coupling constant, a a (1.73GeV) = 0.32^3 (exp) ± 0.016(theo) ||, we get 
0.13GeV ^ fi c ^ 0.31GeV, which is about the same order as the one estimated by the instan- 
ton induced four-Fermi interaction |33|,|35|] or by general effective four-Fermi interactions |32| . 



But, this should be taken as an order of magnitude, since for such a small chemical potential 
the higher order terms in 1/fi expansion, which we have neglected, are as important as the 
leading term. 

So far we have not included the temperature effect in analyzing the gap. The temperature 
effect is quite important to understand the heavy ion collision or the final stage of the 
evolution of giant stars because quarks will have to have high energy to bring together to 
form a super dense matter. The super dense and hot quark matter will go through a phase 
transition as it cools down by emitting weakly interacting particles like neutrinos. 

At finite temperature, T, the gap equation ( |5.35| ) becomes, following the imaginary-time 



formalism developed by Matsubara fl3"6" 



q 2 

Us rp 

9tt , 



E 

l = — oo 



dg 



A(uj r , 



2ttuI + A 2 (w n ) +q- 



In 



A 



(6.6) 



where u) n = %T(2n+l) and q = vp-q. We now use the constant (but temperature-dependent) 
gap approximation, A(u; n ) ~ A(T) for all n. Taking n' — and converting the logarithm 
into integration, we get 
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V ; 18tt J±l 2nJo 



A(T) 



w 2 +A 2 (T)+g 2 x + K-^o) 2 
Using the contour integral (see Fig. 4) |37|, one can in fact sum up over all n to get 



367T 



-T [ dq [ dx (f — 

■ J Jo 2ni Jc 1 



du 



1 



+ e" w / T (cu 2 - g 2 - A 2 (T)) [(u; n - iu ) 2 + x] 



(6.7) 



(6.8) 



Since the gap vanishes at the critical temperature, A(Tc) = 0, after performing the contour 
integration in Eq. (|6.8|) , we get 

(iTT c ) 2 + x-q 2 tanh[g/(2T c )] 



36tt 2 



A 2 

dq I dx 



{ [(ttT c ) 2 + x - g 2 ] 2 + (2vrT c g) 2 2q 
(ttTc) 2 + q 2 -x coth [Vi/(2r c )] 1 



' [(ttT c ) 2 + g 2 - a;] 2 + (2nT c ) 2 x V? j ' (6 ' 9) 

At high density A 3> T^, the second term in the integral in Eq. ( |6.9D is negligible, compared 
to the first term, and integrating over x, we get 



9l 



367T 2 Jo 



Ac tanh?/ 
dy 

y 



In 



A: 



(vr/2) 2 + r 



367T 2 
,2 



1 , tanh-u , , 9 r Xc , tanh-u , A 2 
dy -lnA 2 + / dy In -f + ■ ■ 

2/ Ji y y l 



(6.10) 



= -f^ [(In A c ) 2 + 2A In A c + const. 

367T Z L J 

where we have introduced y = q/(2Tc) and A c = A/(2Tc) and A is given as 

A= [ 1 dy t -^y+ fdy tan ^- 1 =lnffl +7 . 
jo y Ji y Kit/ 

Therefore, we find, taking the Euler-Mascheroni constant 7 ~ 0.577, 

„A 



T. 



c 



—A ~ 1.13 A. 
2 



(6.11) 



(6.12) 



As comparison, we note in the BCS case, which has a contact four-Fermi interaction with 
strength g, the critical temperature is given as 

tanh z 



g I dz 

la z 



dz r 1 , tanh z 
— + / dz 

1 z Jo z 



9 

ji 

g In ( e A Cj 



00 , 1 — tanhz 
dz 



(6.13) 



where cj c (3> 1) is determined by the Debye energy, Cj c = ujr,/{2T c ). It shows that the ratio 
between the critical temperature and the Cooper-pair gap[] in color superconductivity is 



same as the BCS value, e 7 /vr ~ 0.57 PP|Bjl. 



7 In the literature, the BCS gap is defined as twice of the dynamical mass, 2A 
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VII. APPLICATIONS 



It is believed that the core of compact stars like neutron stars may be dense enough 
to form quark matter and may shed some lights on understanding QCD at high density. 
The properties of compact stars can be investigated by studying the emission of weakly 
interacting particles like neutrinos or axions, which is the dominant cooling process of the 
compact stars 0,0. Since the emission rate depends on the couplings of those particles, 
it is important to understand how the interaction and the coupling of neutrinos or axions 
change in dense quark matter. 

Neutrinos interact with quarks by the exchange of neutral currents, which is described, 
at low energy, as four-Fermi interaction, 



C 



vq 



G 



(7.1) 



where G F = 1.166 x 10 5 GeV 2 is the Fermi constant. Again, by decomposing the quark 
fields as in Eq. (p.5|) and integrating out ip- modes, the four-Fermi interaction becomes 



C 



G 



vq 



7^= 2 ^+l(vf, x)ip +L {v F , x)v L {x)fv L {x) + 



(7.2) 



where the ellipsis denotes the higher order terms in the power expansion of l//i. Since the 
four-fermion interaction of quarks with opposite momenta gets enhanced a lot at low energy, 
as we have seen in the previous section (Sec. |V|), it may have significant corrections to the 
couplings of those weakly interacting particles to quarks. We first calculate the one-loop 
correction to the neutrino-quark four-Fermi coupling by the marginal four-quark interaction: 



Gi 



8£vq = -^-ip+ L {vF,x)llj +L {v F ,x)u L (x)fu L (x) 



X 



m_ S A 

2M 2 tv;us 



^(^f, yh^si^F, y)^v{-^ F i vh°^u(-^ v) 



(7.3) 



\%'^ ) +i£ Sf > x ^+l(v F , x)u l (x)V ■ Jv L (x), 



where v F and if F are summed over and g% is the value of the marginal four-quark coupling 
at the screening mass scale M. Now, in fact, because of the kinematical constraint due to 
the presence of the Fermi surface, only the cactus diagrams, shown in Fig. 5, contribute to 
the coupling corrections, which can be summed up as 



5C 



vq 



4G F 
3V2 



E 

V F L 

9?, 



93 



V>+ L (v F , L (v F , x) v\ (x) v L (x) 



Similarly, for axions which couple to quarks as 
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2vr -03, 
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where a is the axion field and fpQ is the axion decay constant, we find the correction to the 
axion-quark coupling to be 



5C 



aq 



3/pq 



E 



V p 



i/j\.(v F ,x)-f 5 ip + (v F ,x) - 



g3V F ■ Va 
2^ + 6-3 



ipl(v F ,x)^ 5 ip+(v F ,x) 



Therefore, as the marginal four-quark coupling approaches 2-rr, the quark-neutrino and 
quark- axion couplings become divergent. Since at low energy g% is quite large in dense quark 
matter, we argue that quark matter will produce neutrinos or axions copiously if the density 
of quark matter is high enough (/x Aqcd) such that the marginal four-quark interaction 
gets enhanced sufficiently 



VIII. CONCLUSION 

We have studied in detail an effective field theory of QCD at high density, constructed by 
integrating out the anti-quarks to describe the low energy dynamics of dense quark matter. 
In the effective theory, the dynamics of quarks is effectively (1+1) -dimensional; the energy 
of quarks does not depend on the perpendicular momentum to the Fermi velocity, which just 
serves to label the degeneracy. At energy lower than the screening mass of electric gluons, 
the relevant interactions are four-Fermi interactions with opposite incoming momenta and 
the coupling with magnetic gluons. 

Because of the dimensional reduction at high density, both four-Fermi interaction and 
magnetic gluon exchange interaction lead to Cooper-pair condensate of quarks in color anti- 
triplet channel for arbitrarily weak coupling. In the ladder approximation, the gap formed 
by the magnetic gluon exchange interaction is much bigger than the one by four-Fermi 
interaction. The color-flavor locking condensate is found to be energetically more preferred 
for high density quark matter with three light flavors, because the color-flavor locking gap is 
bigger than the spin one gap, if we include the four-Fermi interactions. Further, including the 
1/fi corrections and the temperature effects, the critical density and the critical temperature 
for color superconducting phase are calculated by solving the gap equation in the HDL 
approximation. 

Finally, we have calculated the corrections to the quark-neutrino four-Fermi interaction 
and to the quark-axion coupling in dense quark matter due to the marginal four-quark 
interaction. And we found the correction is quite significant and thus dense quark matter 
will copiously produce neutrinos and axions. 
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FIGURES 




FIG. 1. The tree-level matching condition. Wiggly lines denote gluons; solid lines, states near 
the Fermi surface; and double solid line, states in the Dirac sea. 




FIG. 2. The one-loop matching condition for a four-quark amplitude. 




FIG. 3. The 2PI diagram in the leading order HDL approximation. 
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FIG. 4. The contour for integration over the Matsubara frequency. The circles are poles. 
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